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Abstract 



We analyse a simple RPV extension of the MSSM, with bilinear R-parity violation in the soft terms and 
vevs, but not between the terms in the superpotential. The model gives two massive neutrinos, and can fit 
all constraints from neutrino data. We show analytically how to compute the lepton number violating La- 
grangian parameters from neutrino masses and mixing angles. Conversely, we numerically vary the bilinear 
couplings as input parameters to determine the allowed ranges that are consistent with neutrino data. We 
briefly comment on the implications of our bounds for low energy LFV processes. 
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I. INTRODUCTION 



The observed solar |l|J^ and atmospheric neutrino deficits [^|JJ] can be explained by flavour-non-diagonal 
neutrino masses. The atmospheric neutrino anomaly is consistent with Arri^^^ ~ (2 — 5) x 10"^ eV^ (with 
sin^ 20atm > 0.88) The vacuum oscillation interpretation of solar neutrino data requires Am'^oiai- ^ 

10-1° eV^, while the matter enhanced MSW solution prefers the range Ato^^j^^. (10-^° - 10"'') eV^ see 
table ^. In the standard three neutrino framework which offers two independent mass differences, there are 
various mass hierarchies which can explain both solar and atmospheric results. These are summarised in table 
||. The oscillation explanation of the solar neutrino problem and the atmospheric neutrino anomaly requires 
two mixing angles and two hierarchical neutrino mass squared differences, namely Arri^^^ 



« ^<tm- 



Experiment 


Am^ (eV^) 


sin2 261 


tan2 e 


Atmospheric 


(2 - 5) X 10-^ 


0.88 - 1 




MSW-LMA 
MSW-SMA 
MSW-LOW 
Vacuum 
Just-so 


(2 - 70) X 10"" 
(0.4 - 1) X 10"^ 
4 X 10-"' - 2 X 10-^ 
(1 - 6) X 10-1" 
(4 - 10) X 10-12 


0.6 - 1 
10-^ - 10-2 
0.7 - 1 
0.5 - 1 
0.5 - 1 


(2-40) X 10-1 
(1-30) X 10-* 
(1-80) X 10-1 
(1-90) X 10-1 
(3-30) X 10-1 


Chooz 


> 3 X 10-^ 


sine < 0.22 





TABLE L Allowed mass squared differences and mixing angles for MSW-LMA, MSW-SMA, MSW-LOW, Vacuum 
and Just-so stand for MSW large mixing angle, small mixing, and low, and for vacuum and just-so oscillation solutions, 
respectively. See e.g. M. We take the most conservative bounds. 



Spectrum 


Solar 


Atmospheric 


Hierarchical 


Am^2 


Am^3 


Degenerate 


Am23 or Amf2 


Amis 


Pseudo-Dirac 


Am23 or Amf2 


Amis 



TABLE II. Different possible regimes and corresponding mass squared difference. 



We assume that the three neutrino masses m^, i = 1,2,3, satisfy mi < m2 < m^. As we will see, in our 
model we only have two non-zero neutrino masses thus our spectra can only be hierarchical or pseudo-Dirac, 
with Awatin = Amf 3 and Arri^^^^j. = Ami2 , Atojs respectively, see table ||. The 3x3 rotation matrix ||] 
(MNS) which rotates from neutrino flavour (/) to mass (to) eigenstates can be parametrised by three rotations: 
Vfm — R27i{d2'i)Ri'i{0i7i)Ri2{&i2) ■ We neglect the Majorana and Dirac phases, and write 



C12C13 C13S12 Si3 

-C23S12 - C12S13S23 C12C23 - S12S13S23 C13S23 
S23S12 - C12C23S13 -C12S23 - C23S12S13 C13C23 



(1) 



Massive neutrinos can be accomodated in the R-parity (i?p) violating Supersymmetric Standard Model 
|p[-p^, where Rp is defined as (_i)3-B+l+2S^ ^-^^ ^ g g^j.g respectively the baryon number, lepton number 
and spin. Rp is frequently imposed on the Supersymmetric SM to forbid renormalisable baryon and lepton 
number violating interactions. However, the phenomenological bounds on B and/or L violation fl^ can be 
satisfied by imposing i? as a symmetry and allowing the lepton number violating couplings to be large enough 



iThe explanation of the LSND [^] data requires a third neutrino mass squared difference. LSND has not been confirmed 
and will be checked by the future MiniBooNe 



to generate Majorana neutrino masses. We assume in this work that the baryon number violating couphngs 
are absent, and concentrate on the "biHnear" lepton number violating interactions in the superpotential and 
soft terms. We neglect the trilinear R-parity violating (RPV) couplings that can appear in the superpotential. 

Neutrino masses in RPV theories have been studied analytically for many years ||l^,|l3|, and recently 
numerically with more care p^-p6[| . The calculation of the neutrino masses to any given order can be performed 
by calculating the relevant self energy diagrams in the mass insertion approximation for the small RPV mass 
parameters. Alternatively, one can diagonalise the RPV mass matrices of particles which propagate in the 
loops. We use the mass insertion approximation, which allows us to identify analytically the contribution from 
the different parameters and thus use neutrino data to constrain the RPV couplings directly. Our purpose in 
this paper is to consider a simple model in which lepton number violation appears as a misalignment of the soft 
terms with respect to the superpotential. We then calculate the neutrino mass matrix in the mass insertion 
approximation. 

Constraints on RPV parameters from neutrino masses and mixings are significant, as typically they provide 
the stronger bounds on the couplings than other low energy lepton-flavor violating (LFV) processes. In this 
paper we present two approaches in order to constrain the bilinear R-parity violating Lagrangian parameters. 
The first one is an analytic approach in which we obtain expressions for the two bilinear RPV parameters 
appearing in the Lagrangian in terms of neutrino masses and mixing angles. The second approach is a numerical 
one which follows that in ref [ p"7|jT^ . The bilinear parameters are varied as inputs and we use neutrino data 
to constrain the allowed ranges of these parameters. For that, since we only have ranges for the physically 
observed inputs, we constrain our parameter space (6 variables) by using the following inputs from neutrino 
data: 

• Ato^^jjj and sin^ 26latm, 

• Awg^jj and sin^ 20sun, 

• sin^chooz. 

Imposing these constraints implies an indirect restriction on the effective mass TOoff that enters the neutrinoless- 
double beta decay amplitude, i.e., |moff| < '^ifn^ \V^^\. As we will show, we have checked in the different 
analyses we did of that the experimental bound on rUeS is always satisfied. 

The paper is organized as follows, in section II we present the model. Section III presents analytic and 
numerical results. Section IV discusses a model which generates neutrinos masses at tree-level and through 
loops considering only i?p-violating parameter contributions from the misalignment of the /i parameters and the 
vevs. In section V we discuss the implications of our results for low energy lepton flavour violating processes. 
We then present a summary of the paper. 



II. THE BILINEAR RPV MODEL 



In a I^p supersymmctric model, the Higgs and the sleptons can mix (they have the same gauge quantum 
numbers), so the down-type Higgs and sleptons can be assembled in a vector Lj — {Hd, Li) with J : 4..1. With 
this notation, the superpotential for the supersymmctric SM with Rp violation can be written as 

W - fr'HaLj + h-'^^'LjlKE! + X'^'^LjQpDl + hP^H^QpU^ . (2) 

The Rp violating and conserving coupling constants have been assembled into vectors and matrices in Lj space: 
we call the usual fi parameter 1^4, and identify the usual e; = /i^, hl'^ = A^^*^, and = A Lower case 
roman indices i,j,k and p,q are lepton and quark generation indices. We also include possible Rp violating 
couplings among the soft SUSY breaking parameters, which can be written as 

+Al'HuQpU: + A''P'LjQpD''^ + ^A^'''LjLKE^ + h.c. . (3) 

Clearly field redefinitions of the Hd, Li fields correspond to basis changes in Lj space and consequently the 
Lagrangian parameters will be altered. Thus, whenever constraints are placed on the Lagrangian parameters 
the basis in which these are valid must be defined. Alternatively, we can construct basis-independent 



Q 



parametrisations of the couplings and constrain these. In this paper we take the second approach using the 
basis independent parameters 5^,6b, constructed in refs. 21|, which are defined in table Il|. 
It is well known that this model can generate AL = 2 neutrino masses from tree-level and loop diagrams. 
There will be a non-zero neutrino mass at tree-level if /i/ ^ (a/^j /ij/ Vj)vi — that is if (/X4, /is, /i2, /ii) 
is not parallel with (t'4, Us, i^i) In the Vi = {vi) = basis, a neutrino v-^ acquires mass 7713 at tree level 
via a "seesaw", with neutralinos playing the role of the heavy Majorana fermion, and the mass fii^ih^ in the 
place of the "Dirac" mass, see figure |l|. 




FIG. 1. Tree-level neutrino mass in the mass insertion approximation. 

In the mass insertion approximation for l^p masses there are three types of loops that can contribute 
masses to the two neutrinos that are massless at tree level. There are the well-known loops involving trilinear 
I^p couplings A or A' at both vertices. Then there is the Grossman-Haber diagram [0, figure ^, with gauge 
couplings at the vertices and I^p provided by mass insertions of the soft masses which mix the Higgses and the 
sleptons. Finally, there is figure ^, which has one gauge coupling and one trilinear/yukawa, and at least one 
unit of lepton number violation provided by a mass insertion. The last two types of loops have frequently 
been neglected in analytic estimates of neutrino masses. In a previous work pG,pl|, two of us showed that 
the analytic estimates in models with lepton number violating masses (= bilinear ) suffer from two linked 
problems: confusion about the interaction eigenstate basis choice in the Lagrangian, and an incomplete set of 
one loop diagrams contributing to neutrino masses. "Basis independent" formulae for all the diagrams can be 
found in ||2^. For a discussion of the construction of basis-independent invariants, see, e.g. [| 20 - 2^] . 



R-parity violating models can be divided into classes based on which of the diagrams make significant 
contributions to the neutrino mass matrix [m,y]ij [ pO| , pl| : in the first case, the loop diagrams with bilinear Ijtp 
are much smaller than the trilinear diagrams. In this case, [mi,]ij is due to the tree mass and the canonical A 
and A' diagrams, and can have a large variety of patterns depending on the MSSM and I^p inputs. This case 
has recently been carefully studied in refs. In the second case, the loops with bihnear I^p due to soft 

masses make significant contributions to [m,y],y, but the bilinears contributing at tree level((5^) are negligible 
in the loops. This is the case we study in this section. To isolate the effect of the soft , we assume that the 
trilinears are negligible. The third possibility is that loops with lepton number violation from Sp_ are important. 
Such loops can contribute [m^]ij ^ {IQ-Kmsusv), whereas the tree contribution is m^'"'^'^ ~ ^^^rasuSY- 

This generates a very large hierarchy among neutrino masses, so this case is unlikely to be realised. We discuss 



this further in section IV 



We take in this section a model where ^ 0, and where 5\,5y are negligible. Such a model could 

arise if the I^p is originated in the soft terms. It is interesting to consider because it shows the contribution 
of S-'g to [TOi,]y , and allows us to set bounds on the (5^. We want to determine which values of ((^Ij, '^b) 
and ((5^, 5^, reproduce a neutrino mass matrix consistent with the data. We do this in two ways. First we 

assume that the neutrino masses and mixing angles are all known, and analytically solve for 5^ and (5b as a 
function of the masses and mixing angles. This is tractable because the (5^, 5b model only has two non-zero 
neutrino masses, so the algebra reduces to two dimensions. 

Secondly, we start with (5^ and 5b as inputs, varied over "sensible" ranges, and numerically determine 
whether the resulting neutrino mass matrix is consistent with observations using the inputs given above. 

There are three relevant diagrams contributing to [mi,]ij in this model: the tree diagram of figure 0, and 
the Grossman-Haber diagram of figure || with I^p at points IV and VI, and with ^p at points I and VI. Exact 
formulae for these diagrams can be found in |2l|| . Our simplifying assumptions generate a neutrino mass matrix 
of the form 

[^uh - ^^TKK + ^'-^UK^B + + mf5),P^ (4) 

where 



\p\V2mi 




ji — B-V-v 
^ ~ \B\V2mi 


\B\ 






ijijk v-\'\''X^-v 





TABLE III. The basis-independent invariants used to parametrise the bilinear 1^ relevant for neutrino masses, 
together with their value in the (z>i) — basis. They are zero if Rp is conserved. (Note that these invariants have signs: 
for arbitrary vectors a and h, a ■ X -h — —h ■ ■ a). 



h, H, -A 
IV VI 



. VIII 
_i m v-i 



FIG. 2. The Grossman-Haber loop. The blobs indicate possible positions for mass insertions. The misalignment 
between fl and v allows a mass insertion on the lepton/higgsino lines (at points I, or VIII). The soft I^p masses appear 
as mass insertions at positions VI and IV on the scalar line. 



m(,o) = \S^\^msusY (5) 



m 



^ IGtt 

Yg^, and Sb are unit vectors in space: 



- —\SB\'msuSY (7) 



S,= iS;,S^^,S;) . (8) 



We have set all the unknown sparticle masses equal to msuSY = 100 GeV, and neglected the mixing angles 
among MSSM particles. This is phenomenologically reasonable since superpartners have not been detected. 
Using the correct dependence on the MSSM parameters would have two effects: m^^^ would no longer be 

related to m^°^ and m^^ as in eqn (||), and there could be additional (probably small) contributions to [m^]y 
with more complex index structure induced by spartner mass differences. 
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FIG. 3. Schematic representation of tiic charged loops with one gauge and a Yukawa couphng. These diagrams occur 
if gauginos mix with charged leptons — that is if S^j, ^ 0. The blobs indicate possible positions for ^ interactions; 
each diagram contains two blobs. The I^p interactions can be trilinears (at VII) or mass insertions. The misalignment 
between pL and v allows a mass insertion on the lepton/higgsino lines (at points I, III, or VIII), or on the scalar line at 
point V . The soft I^p masses can appear as a mass insertion at position IV. 



III. RESULTS 
A. dfj, and 5b bs a. function of masses and mixing angles 

In this section, we wish to go from the data to the Lagrangian parameters 5^ and 5b- We assume that 

the three neutrino masses (one of which is zero) and three mixing angles are all known exactly; from these 
six inputs we wish to analytically solve for the six parameters (5^ and 5^, (i = 1,3). The data do not fully 

determine 5^ and Sb, because these two vectors are not required to be orthogonal. An orthonormal basis would 
be 5fj, and 5b j_ where 

--{5% ^,5^,^,5], J , (9) 



and 

s'b^ = - i^B ■ s^)s; . (10) 

The parameters 5^ and 5b can be expressed in terms of the two neutrino masses and the mixing angles, and 
one additional parameter which could be taken as the angle between (5^ and 6b (which we will call p). 

We first express the neutrino mass matrix [m^]ij in terms of the orthonormal basis {5^, 5bj^)- Since [m,y]ij 
is a mass matrix in the two dimensional space spanned by 5^ and fe, it has two non-zero eigenvalues and the 
problem rc;(luccs to two dimensional algebra. Secondly, we diagonalise [m^]ij, and compare the eigenvalues and 
eigenvectors to the observed masses and mixing angles. This allows us to determine 6^ and (5^ as a function of 
the masses and mixing angles, and one free parameter. The free parameter arises because (5^ and 6b are not 
orthogonal, so what fraction of the heavier neutrino mass is due to is a free parameter. 

In the ((5^, 5b±) basis, can be written 

[m^]ij = (m^°) + 2m^°l cosp + m^^ cos^ p)6l6j^ 

+ {m^^l sin p + m^^hmp cos p){5id^s^ + d%Jj^) + m^shm'^ p6i^J's^ (11) 
= m^6'Jj, + m^Bi6'J'B^ + S^bJI) + ms^jj^^ij^ (12) 



where p is the angle between 5^ and 5b '■ cosp = 5^- 5b- The matrix [m^]ij has two non-zero eigenvalues, so 
the neutrino mass spectrum is either hierarchical with Am^^^ — m\ and Am^^^ — m|, or pseudo-Dirac, with 
Am^j„ = m§, and Am^^; = 7713 — m|. The eigenvectors are the last two columns of V/,„ in eq. (|l|), and lie 
in the 2-d sub-space spanned by (J^j, (5^, (5^) and (5^, 5'j^,5^^). 

We can now determine (<5]j, 5^g,5'^) and (5^, (S^, (5^) as a function of the neutrino masses and mixing angles, 

and the angle 7 between (^,(5^,(5^) and (V13, V23, ^33)- We solve for the angle p (between 5b and 5^) as a 
function of 7 in equation (17). 

In the mass eigenstate basis spanned by V/2 and V/3 (/ : 1..3), the neutrino mass matrix is 



TO2 

TO3 



(13) 



Suppose that ((5^, 5'^^, 5'^^) is rotated by an angle 7 with respect to (V13, V23, V33), so that the mass matrix (|l 
is 





TUB 


m^B 




m^B 





cos^ 7TO2 + sin^ 77713 cos 7 sin 7(7772 — 7^3) 
cos 7 sin 7(7772 — 7773) cos^ 77773 -|- sin^ 77772 

in the basis of {5jj_, 5f^,5f^) and {5g^ , 5g^ , 5^^). By comparing eqn ( |l^ ) to eqn (|l^), we find 

(0) .9 ^iri2 -2 / 2 -2 \ 

777^' sm p=:|-^|dB| 77J5(7sysm p = (cos 77712 + sm 77713) 



,(0) 



cot p 



16n' 



cos 7 sin 7(7712 — m^i) 



fl-— ) 
V 167r/ 



77127713 



cos^ 7?772 + sin^ 77773 



7772 7773 



cos^ 77772 -l- sin^ 77773 y 167r (cos^ 77772 -I- sin^ 77773)^ 



(14) 

(15) 
(16) 

(17) 



T7-^5l^ cosjVfs + sin7V/2 (18) 
1 , 

|7 — I ^B^ = - sin 7^/3 + cos 7I92 ■ (19) 

We can now determine the values of {5j^,5f^,5^) and {5B,5g,5g) corresponding to atmospheric and different 
solar solutions. 

We first compute values of and \5b\ that give solar and atmospheric neutrino masses. 
To get degenerate masses 7772 ~ 7773 ^ \J Atti^^^, the angle p between 5b and 5^l must be of order ±7r/2. So 
5b must be approximately orthogonal to 5^, and ~ |i5^p ~ ^ l\rr?^^^lmsusY ■ This seems unnatural 

and fine-tuned, so we do not pursue this possibility. 

For the case of hierarchical neutrino masses, we plot in figure^the allowed values of and \5b\ corresponding 
to combined constraints from Chooz and atmospheric (SuperK) results and three of the solar solutions. The 
allowed parameter space corresponds to the space contained between the lines in figure |4| The shortest pair 
of lines corresponds to the LMA MSW solution, the middle lines to the SMA solution and the longest pair of 
lines to the Vacuum solution. We obtain these two lines by varying the angle 7 for the maximum and minimum 
allowed values of the masses m2 and 7773, as given in table |^. In this plot (fig.p, the left [bottom] end of the 
allowed parameter space corresponds to 7773 induced by (5^ [5^]. 

It was pointed out in ||l^ that their MSUGRA bilinear model was effectively controlled by two vectors, 
5y, and the misalignment between p and the trilinears (= 5\). The parametrisation of equations 
should therefore approximately describe the results of p6| , replacing 5b \yy 5\. Their heavy neutrino mass 7773 
is induced by the tree diagram of figure ^ so the angle 7 (between the eigenvector corresponding to 7773 and 

is small. It is easy to see from equation (^8|) that the atmospheric mixing angle is therefore controlled by 
the relative size of 5^ and (5^, as was indeed found by We see from equation (|l^) that the solar mixing 
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FIG. 4. Allowed region in parameter space for various solar oscillation solutions combined with the atmospheric and 
Chooz constraints. The allowed region is between the lines; the longest pair of lines correspond to the vacuum solar 
solution, and the shortest [middle length] lines to the LMA [SMA] MSW solution. 



angles are determined from the relative size of the elements of Sb^, rather than from the elements oi Sb- This 
means that the solar angle will sometimes show a correlation with the elements of (5b, specifically if 6b is 
approximately perpendicular to S^. This agrees with the results of | |T^ . 

The Chooz result^ require Ve3 = cosjS'^ — sin 7(5^^ ^ .1. For large angle solar solutions (LMA), Ve2 ~ 
l/-\/2, so 5^ ~ sin 7, and 6'^^ ~ cos 7, and Ves is small due to a cancellation between 6^ and 6'^^. For 

small angle solar solutions, both 6^ and 6%^ may be small. This scenario has cosmological attractions, because 
it would allow a baryon asymmetry produced before the electroweak phase transition to survive. Baryon + 
lepton number violating processes are in thermal equilibrium before the electroweak phase transition, so a 
primordial asymmetry must be stored in a conserved quantum number (which could be one of the i?/3 — Li) 
to survive. For B/3~ Le to be effectively conserved before the phase transition, we need |^ {B^n'^ — n^B'^) < 
2 X 10-^|M||5|/tan/5. | For 6^ ~ 10 ®, Sb ~ 10 ^, as could generate the atmospheric and SMA solar masses, 
this is satisfied for (5;;; < .1 or (5| < 2 x IQ-^. 

In the case of a hierarchical spectrum, then only for very small values of the angle 7 do we find that 
cos^ 7m2 sin^ 77713. In this case \5^\ will generate the neutrino with a mass — > -^/Amf^ and \Sb\ 
will generate the neutrino with a mass m2 \/ Am^^^^^. For all other values of 7 then |(5_b| \J ^"^atm ^"^^ 



^The original two-neutrino analysis of the Chooz data combined with the SuperK best fit atmospheric value gives 
|t/e3|'^ ^ 4. 10"'^. In the framework of three-neutrino mixing, the analysis of Chooz combined with the SuperK atmo- 
spheric data gives [f/esi' ^ 2. lo-^ see ref. |Q and references therein. 

''We have added a factor of tan/3 with respect to the first reference of |2^. This comes from requiring that the lepton 
number violation due to the soft masses be out of equilibrium; the minimisation conditions of the scalar potential say 
that the RPV soft masses vn?Ai = — -Bj tan/3, so can be larger than the B terms for large tan/3. 



^ a/ Am'j^i^^. That is, the loop is giving the largest contribution to the neutrino masses. So the "measure" 
on parameter space determined by the angle 7 does not correspond to natural theoretical expectations. For 
large values of 7 the atmospheric mixing angle is now controlled by the relative size of Sg^ and S"^^ . 

B. Masses and mixing angles from and Sb 

In this subsection we present numerical results of our scan of parameter space varying the inputs pa- 
rameters in the range of |(5^| < 10^^ and |(5_b| < 4.10^^. This scan allows us to consider the case when tree 
contributions are larger, of the same order or smaller than the loop contributions. We present in figure ^ 
the allowed region in parameter space in the \6b\ vs |(5^| plane, these results are in agreement with the ones 
obtained in ref. pl[ |. The different points displayed correspond to combined constraints from Chooz+SuperK 
atmospheric data and one of the solar solutions. As expected, figures ^ and ^ are quite similar. 

Solan- CH 00Z+ SuperK 

1 1 1 1 1 1 1 1 1 1 
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I I I I I I I I I I I 

LO"' 10"' 

^ 

FIG. 5. We present the available region in parameter space for a combined fit using Chooz, atmospheric and one of 
the solar oscillations solutions. 



In order to study the spectrum for each of the possible combined solutions we give samples of the different 
mass spectra for different values of the input parameters (5^ and 6b- For illustration, in figures ^,0, || and ^ we 
present the obtained spectra for the combined constraint of Chooz and the atmospheric SuperK data with one 
of the solar solutions as functions of \6^\ and \6b\ ■ We also plot m-e// and rrii and we see that the bounds 
from neutrinoless double beta decay and the cosmological constraint rrii ^few eV are respectively 
fulfilled. From these plots one can read that the scan shows us only hierarchical spectra for all combined 
constraints. 
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FIG. 6. The mass spectrum is represented together with rrias and the sum of the eigenvalues in the case of a combined 
fit with Chooz, Vacuum solar solution and atmospheric SuperK constraints. 



2| and [T^ we plot the difFerent solar and atmospheric mixing angles as functions of the 
R-parity violating parameters -jj-r ^-nd jg-^- The allowed ranges for the RPV parameters agree with the 

estimates of equations and ( jl^ ) . 

In general, the observed discreteness in all our plots is in part due to the step taken in our scan of parameter 
space and we note that the dependence on a single parameter is complicated since many solutions can be found 
for the same input value. Also the same value of the physical parameter involved can correspond to difFerent 
values of the inputs. 
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IV. MODEL WITH ONLY S^, 

In this section, we show that a realistic neutrino mass matrix cannot be generated in an RPV model 
where only the invariant 5^ is non-zero. This toy model is interesting because it allows us to study whether 
the 0(|(5^p) loop contributions to [m^]ij need to be taken into account. We assume that the only non-zero 
invariant is S^, and ask whether we can generate m2 from the (5^ loops of figure ^ if is due to the tree level 
diagram of figure 0. The answer is no, even for large tan/3, because the loop mass is necessarily small (of order 
the vacuum or just-so masses), and the solar mixing angle is also small. 

The analysis of this section is quite relevant as this type of model has been constrained using other low 
energy LFV processes, as in | |2^ . However, as mentioned in the introduction, typically neutrino data will put 
tighter constraints on the RPV parameters and as we have seen it is not possible to accomodate all neutrino 
data simultaneously. 

We take 

\S^\ ^ (20) 
msuSY 

and 3> so that the heaviest neutrino mass TO3 explains the atmospheric neutrino anomaly and is 
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FIG. 7. The mass spectrum together with the effective mass rrics and the sum of the eigenvalues are represented for 
a combined fit with Chooz, Superk and MSW-LOW oscillation solutions. 



consistent with Chooz. The solar neutrino mass will come from the "see-saw-type" loops of figure ^ with one 
mass insertion on an external leg. 
The diagrams of figure ^ can have R-parity violating mass insertions at (I,V), (1,111) or (III, VIII). The 
potentially largest is (I,V) (see for exact formulae), because it is proportional to tan/3: 

2 

TTl 

(/, V) ^ gS^Sj, tan f3-r^ . (21) 

We refer to this as "see-saw-type" because it is a mass between the neutrino who acquired mass at tree 
level, and the massless-at-tree- level neutrinos 1^2 and vi . Neglecting the muon and electron masses with respect 
to the tau mass, equation ( pi] ) is a mass between and Ut-. The largest see-saw mass it could induce for 1/2 
would arise if 1^2 ^ i^r, and is of order 

^^^^,_/tan^y . (22) 



This could generate m2 ^ lO^^m^J''^'^ for tan/3 50. So if all free parameters are stretched, it might just be 
possible to get a "see-saw" loop mass of order the Vacuum or JUST SO solar neutrino masses. 

The neutrino which gets a loop mass will only have a very small Ve component, so will not have the 
large mixing angle required for the Vacuum and JUST SO solar solutions. However, recall that we choose 
SJ^ ^ 5^ ^ 6'^, so that the tree level mass explains the atmospheric neutrino anomaly and agrees with Chooz. 
The loops mix 1^3 with so the second massive neutrino will be composed largely of v-r and v^, and the 
massless neutrino will be mostly Ve- An Rp violating model with 6fj, ^ and all the other invariants zero 
therefore can not explain observed neutrino data. 
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FIG. 8. The mass spectrum together with the effective mass rrics and the sum of the eigenvalues are represented in 
the case of a combined fit with Chooz, Superk atmospheric and MSW-SMA oscillation solutions. 



V. EFFECTS ON LOW ENERGY LFV PROCESSES 

It is also of interest to comment on the effects of the allowed values for the different dg , Sj^ from neutrino 
masses and mixings on low energy lepton flavour violating processes. In this model there can be tree level 
contributions from the RPV parameters to some of the AL = processes. We can estimate the contributions 
to processes like fi —^ ej, fi — > eee /i — e conversion in nuclei, and other LFV processes at low energy. 

In general the branching ratios for these processes will be proportional to 



(23) 



Thus it is trivial to see using simple estimates that the branching ratios of low energy lepton flavour 
violating processes will automatically be suppressed given the constraints from neutrino data. The obtained 
branching ratios are well below current and anticipated experimental values. 

Our simple model for generating neutrino masses and mixings correlates the ratios of branching ratios 
of different LFV low energy processes. The exact values of these correlations will depend on the combined 
constraint of solar, atmospheric and Chooz experiments which is used. 

A generic consequence of this model is that the lightest supersymmetric particle (LSP) will decay. The 
details of this analysis will be presented elsewhere. 



VI. SUMMARY 



In this paper, we studied a minimal model of bilinear Rp violation which is consistent with the atmospheric 
and solar neutrino anomalies. In the "basis-independent" formalism in which we work, the model has (5^ ^ 
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FIG. 9. The mass spectrum together with the effective mass rricft and the sum of the eigenvalues are represented for 
a combined fit with Chooz, Superk atmospheric and MSW-LMA constraints. 
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FIG. 10. The allowed atmospheric and solar angles for a combined fit of MSW-LMA, Chooz and SuperK atmospheric 
data. The atmospheric (solar) angle is plotted as a function of |i5j| and \Sf^\ (j^Sj^l and in (a) and (b) (in (c) and 

(d)) respectively. 



6b 7^ 0, and all other invariants zero In the basis where there are no ^.p trilinears, this model has R-parity 



*The vectors and lower case indices are in leoton flavoun aoace: our notation is exolained in the introduction and in 
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FIG. 11. The allowed atmospheric and solar angles for a combined fit of MSW-SMA, Chooz and SuperK atmospheric 
data. The atmospheric (solar) angle is plotted as a function of |i5^1 and \S'^\ {\Sg \ and \5% |) in (a) and (b) (in (c) and 
(d)) respectively. 



violating BiHuL^ masses and sneutrino vevs. Such a model could arise if the superpotential was Rp conserving, 
but R-parity violation was introduced through the soft terms. 

In the limit where we neglect MSSM mixing angles and set all spartner masses to a common scale msuSYi 
the neutrino mass matrix is controlled by the 6 parameters (5^, (i: 1..3). There are two massive neutrinos. If 
we make the simplifying assumption that we know the neutrino masses and the mixing matrix (MNS) exactly, 
we can solve for 5^ and 5b as a function of one additional parameter, which can be taken as the angle between 
5a and 5b- We plot the allowed values of jJ^I versus \5b\ for different solar solutions in figure ^. See also figure 
|5|. This gives some idea of the weak-scale I^p Lagrangian mass terms which can induce a phenomenologically 
correct neutrino sector. In reality however, we do not know the neutrino masses or the MNS matrix elements 
exactly, so we numerically vary the six inputs parameters, and find ranges where the neutrino masses and MNS 
matrix elements are consistent with combined fits of SuperK, Chooz and solar data. 

This model differs from usual models of bilinear I^p in that it has different loop diagrams. 5b generates 
a loop contribution to the neutrino mass matrix of order g'^\5B\'^'msuSY via the Grossman-Haber diagram 
of figure H This diagram has frequently been overlooked in the literature. If (5^ cx 5b, then this is a loop 
correction to the tree mass and so its not useful to generate another neutrino mass. However we treat 5p, 
and 5b as phenomenological weak-scale inputs, and allow them to vary independently. The Grossman-Haber 
diagram is therefore the source of our loop neutrino masses. This is in contrast to a MSUGRA bilinear model, 
where /x is misaligned at the GUT scale with respect to the trilinears, and the loop diagrams are proportional 
to trilinear couplings squared. The trilinear diagrams are more important than the Grossman-Haber diagram 
in MSUGRA models, because renormalisation group running makes the latter a loop correction to the tree 
level mass [5^ ^ 5b)- 

The simplicity of this 5^ — 5b model is attractive: neutrino oscillations can be explained with six additional 
parameters (and all of the MSSM!). This is in contrast to most models with trilinear R-parity violation, 
where a multitude of couplings contribute to the neutrino masses. The phenomenology of this model is quite 
different from other RPV violating models with respect to the scalar sector of the theory. For neutrino masses 
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FIG. 12. The allowed atmospheric and solar angles for a combined fit of MSW-LOW, Chooz and SuperK atmospheric 
constraints. The atmospheric (solar) angle is plotted as a function of and \S'^\ {\Sg \ and \S% |) in (a) and (b) (in 
(c) and (d)) respectively. 



y/ Am'^^i and a given MNS matrix, the mass terms Bi and /ii in the Lagrangian 



7713 = v^rn^tm and "^2 . _ 

are determined as a function of one free parameter, which is the angle 7 in equations (|l^) to (|19|). This 
parametrisation can describe more general models of bihnear under certain circumstances, as was discussed 
in section [II A. We have also shown that combined solutions of SuperK, Chooz and all possible solar oscillations 



solutions can be obtained. 
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